Reconstruction of supernova z/„, u T , v„, and v T neutrino spectra at scintillator detectors 
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We present a new technique to directly reconstruct the spectra of v T , v^, and v T from a 
supernova, using neutrino-proton elastic scattering events (y + p — > v + p) at scintillator detectors. 
These neutrinos, unlike v e and v £ , have only neutral current interactions, which makes it very 
challenging, with any reaction, to detect them and measure their energies. With updated inputs 
from theory and experiments, we show that this channel provides a robust and sensitive measure of 
their spectra. Given the low yields and lack of spectral information in other neutral current channels, 
this is perhaps the only realistic way to extract such information. This will be indispensable for 
understanding flavor oscillations of SN neutrinos, as it is likely to be impossible to disentangle 
neutrino mixing from astrophysical uncertainties in a SN without adequate spectral coverage of all 
flavors. We emphasize that scintillator detectors, e.g., Borexino, KamLAND, and SNO+, have the 
capability to observe these events, but they must be adequately prepared with a trigger for a burst 
of low-energy events. We also highlight the capabilities of a larger detector like LENA. 
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I. INTRODUCTION 

The detection of neutrinos from a core-collapse super- 
nova (SN) is the key to understanding SNe and neutrino 
properties. The fact that SNe emit all six flavors of 
neutrinos and antincutrinos, with average energies that 
depend on their different cross sections, offers a rich 
potential to reveal the detailed properties of the proto- 
neutron star. And the fact that these extreme conditions 
can lead to strong neutrino mixing may amplify the 
effects of neutrino properties too subtle to be seen in the 
laboratory. However, to fully separate and identify the 
astrophysical and particle physics effects, all the flavors 
must be detected and their spectra measured. While 
SN 1987A provided evidence of the expected strong 
neutrino emission from a core-collapse supernova [l|, Q , 
only ~ 20 v e events were detected, allowing only modest 
precision in the spectrum measurement and no clear 
evidence of what mixture of the initial v e and v^/y T led 
to the received v e 

It will be challenging to experimentally measure and 
theoretically interpret supernova neutrino data. The 
prospects for success depend crucially on the values of 
the emission parameters and especially on the differences 
between flavors, e.g., their average energies. Compared 
to a decade ago, it is now widely thought that the 
average energies and the differences between flavors 
are both less IZHI (see [H for reviews), which 
reduces the expected numbers of events and the effects of 
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neutrino mixing, making the challenges greater. Further, 
there is now an appreciation that neutrino mixing 
is much more complicated due to neutrino-neutrino 
interactions p~2l - [25| (see [lH for a review), making the 
importance of high-statistics measurements with flavor 
and energy dependence even greater. 

It is likely that a Milky Way supernova will occur 
in the coming decades, and it is exceedingly important 
that we are prepared to capture the most detailed data 
possible. Present detection capabilities are much better 
than for SN 1987A, and the likely supernova distance 
much closer, so a large detected yield is expected. The 
v e spectrum will be measured in Super-Kamiokande and 
other detectors with ~ 10 4 events. The v e spectra could 
be measured well in one of the proposed liquid Argon 
detectors. For these flavors, there are charged current 
detection channels with large cross sections and good 
spectral fidelity. 

The other flavors, v^, v T , v^, and v T , collectively called 
v x and assumed to be similar, while harder to measure, 
are especially important because they constitute the 
bulk of the emission and drive neutrino mixing effects. 
Because the charged current channels are energetically 
forbidden, detection depends on the smaller neutral 
current cross sections; this is partially compensated by 
the correspondingly higher emitted average energy and 
the four flavors. The most serious difficulty is separating 
these events from other channels and measuring their 
spectrum without having charged leptons in the final 
states. 

The larger the average energies of all flavors and their 
spectral differences, the easier it is to measure v x spectra. 
Ideally, the charged current channels would show two 
distinct spectral components, due to mixing. However, 
on the basis of SN 1987A data and contemporary 
supernova simulations, this seems unlikely. Some neutral 
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current channels, e.g., those on 12 C and 16 0, lead to the 
emission of gamma rays with energies characteristic of 
the nuclei and not the incident neutrinos [13, • The 
yields are low for what now appear to be reasonable 
average energies, and sensitivity to the assumed average 
energy is degenerate with uncertainties in the total flux 
and the nuclear cross sections. Neutral current neutrino- 
electron elastic scattering does yield spectral information, 
but it is very difficult to isolate these events from charged 
current events in this and other channels. 

As a solution to these problems, Beacom, Farr, and 
Vogel (BFV) pointed out that neutrino-proton elastic 
scattering (y + p — > v + p) in scintillator detectors 
could give a large yield of separable v x events with 
spectral information, provided that the low-energy de- 
tector backgrounds are low enough [29j . Fairly opti- 
mistic assumptions about the properties of then-future 
detectors and the v x average energy were required. This 
detection reaction mainly probes high neutrino energies, 
so it is important to ask if it is viable with the known 
capabilities of present detectors and the lower average 
energies assumed today. 

We show that this technique is indeed viable with 
realistic inputs and in fact is essential to adequately 
understand supernova emission and neutrino properties, 
especially in light of collective oscillation effects [2r| . 
We provide new detailed calculations with contemporary 
inputs for several detectors: the presently-running Borex- 
ino and KamLAND, the near-term SNO+, and the much 
larger proposed LENA. These results are needed so that 
the experiments have triggers in place that will ensure 
that this data is not missed, and to show how to interpret 
the signals. Most importantly, we develop a new method 
to directly invert the measured proton spectrum for the 
unknown v x spectrum, allowing one to go beyond the 
thermal spectra assumed by BFV. 

The outline for this paper is as follows. In Sec. [H] 
we discuss why learning about SN neutrino spectra is 
a complicated problem, and how measuring v-p elastic 
scattering events addresses the issue. In Sec. IIII1 we 
review the general framework of v-p elastic scattering at 
scintillators. In Sec. IIV1 we present the expected signals, 
and, in Sec. El our prescription to reconstruct the v x 
spectrum at Earth. We discuss, in Sec. I VII possible 
improvements if a larger detector is built. We discuss 
phenomenological implications in Sec. I VIII and conclude 
in Sec. IVTiTl 



II. COLLECTIVE NEUTRINO MIXING 

Ideally, one would like to determine the SN neutrino 
emission parameters, study SN explosion dynamics, and 
determine neutrino parameters from a measurement of 
SN neutrino spectra. Different aspects of the uncertain 
SN astrophysics and neutrino oscillations are coupled due 
to the effects of neutrino mixing in dense matter. These 
effects can be strong, have a complicated phenomenology, 



and are still not completely understood. All this 
makes the problem of reconstructing the parameters of 
SN emission and neutrino mixing a highly challenging 
problem. 

While neutrino parameters may eventually be de- 
termined from cosmology [30[ or oscillation experi- 
ments 13 lj . there is no other way to probe SN neutrino 



emission or explosion parameters in detail. One needs 
measurements of the received neutrino fluxes, which, 
when interpreted with guidance from state-of-the-art SN 
simulations [HI [H| , could reveal important aspects of SN 
physics. 

This inverse SN neutrino problem was already difficult 
in simple neutrino mixing scenarios |32| . However, 
there was at least one major simplification. The 
flavor evolution did not depend on the initial spectra 
themselves, i.e., the Hamiltonian was independent of 
the neutrino fluxes. As a consequence, observing the 
v x spectra wasn't necessary for determining the flavor 
evolution of v e or v e . 

However, it has since been realized that SN neutrinos 
are also subject to so-called "collective effects," due to 
neutrino-neutrino forward scattering [l2l - [25j . This is 
unavoidable near the SN core, and can lead to large 
effects for almost the entire duration of the SN burst. 
In the presence of these collective effects, the flavor 
histories of all species get coupled to each other along 
their trajectory. The evolution depends explicitly on the 
flavor-dependent fluxes of the neutrinos. Therefore, to 
calculate the flavor evolution of any neutrino flavor in 
a SN, one requires knowledge of initial conditions for 
all others - a situation fundamentally different from the 
older paradigm. However, if one knows the final state 
of all species, one can choose possible initial conditions 
for the neutrinos, evolve them forward, and check if they 
reproduce the final spectra. This allows one to interpret 
SN neutrinos. Of course, this bootstrap is meaningful 
only if one has a complete characterization of the final 
state. Without that, there will be strong degeneracies - 
a large suite of initial conditions can produce the same 
final spectra for appropriate choice of neutrino and SN 
emission parameters, and no firm inference is possible. 

In this light, it becomes crucial that all SN neutrino fla- 
vors be observed at Earth. This can be achieved through 
a detection of v-p elastic scattering events in addition to 
the charged current measurements. Detecting all flavors 
would break degeneracies and allow determination of 
the primary neutrino spectra. This would also allow 
robust analysis of model independent signatures of SN 
neutrinos, e.g., to determine neutrino mixing parameters 
or probe Shockwave dynamics [33l . |34| . It is with this 
motivation that we study the detection of v-p clastic 
scattering events from SN neutrinos. 
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III. GENERAL FRAMEWORK 

A core-collapse SN provides a neutrino fluence, i.e., 
time-integrated flux, dF/dE at Earth, spread over an 
energy range E ~ (5 — 50) McV. These neutrinos inter- 
act with N p free protons in the detector through neutral 
current elastic scattering. The differential cross section 
da/dT is broad and slightly forward peaked, leading to 
proton recoils with kinetic energies T < 5 MeV. Each 
recoiling proton, as it is brought to rest, deposits energy 
in the scintillator. These protons are slow, so the 
scintillated light is quenched, i.e., the number of photons 
scintillated corresponds to a lower effective proton recoil 
T' < 2 MeV. One observes the effective proton spectrum 
dN/dT', and the objective is to extract the neutrino 
fluence dF/dE. 

The observed event spectrum due to SN neutrinos is 



dN N p 
dT> ~ dT'/dT J E 



. „ dF da , „. 
dE- — —(E) 
dEdT y ' 



(1) 



A neutrino of energy E can produce a proton recoil 



energy between and T n 



2E 2 /rn p , where m p is 



the proton mass. In other words, a minimum neutrino 
energy E m - m = -J m p T/2 is needed to produce a recoil 
energy T. The recoiling proton is unbound from its atom 
and molecule, so that its energy loss in the medium is 
dominated by collisions with electrons. 

There are three ingredients needed to calculate the 
time-integrated SN neutrino event spectrum due to 
v-p elastic scattering at a scintillator detector: (i) 
the neutrino fluence dF/dE over the signal duration 
At, (ii) the cross section da/dT, and (in) detector 
specific information, e.g., the number of target protons 
N p , quenching function T'(T), and energy resolution. 
Additionally, the relevant backgrounds over the duration 
of the burst are needed to estimate signal significance. 



A. SN neutrino fluence 

A SN emits a total energy £ w 3 x 10 53 erg over 
a burst of A( « 10 s in neutrinos of all six flavors. 
The neutrino flavors v^, v T , and their antiparticles, have 
similar interactions and thus similar average energies 
and fluences. Therefore, the total energy is divided as 
£ = £ Ve + £jj e + 4£„ x . The various flavors have different 
average energies - lowest for v e and highest for v x . The 
fluence in each flavor is distributed in energy according 
to a normalized spectrum dip a /dE. A SN at a distance 
d from Earth thus provides a neutrino fluence 



dF 
dE 



E 



i 



dF n 



dE 4ird 2 



E 



£g dip a 
(E a ) dE 



(2) 



where dF a /dE is the fluence in each flavor. For a 
neutral current process, the initial flavor distribution 
is not important: only the sum of all active-flavor 



neutrinos is relevant. Assuming no active-sterile mixing, 
one can ignore subsequent neutrino oscillation effects for 
calculating this proton recoil signal. 

The parameters of the neutrino fluences are not 
known accurately, and the objective is to measure them. 
For concreteness, we choose as a nominal spectrum 
d(p a /dE = (128/3) (E 3 /(E a ) 4 ) exp(-4E/{E a )). This is 
one variant of the Keil parametrization for the spec- 
trum [35|], and is normalized as J dE d<p a /dE = 1. The 
fluence in each flavor is then 



2.35 x 10 13 £ a E 3 



dF, 



dE cm 2 MeV d 2 (E a ) 



exp 



AE 



(3) 



In the last expression, £ a is in 10 52 erg, d is in 10 kpc, 
and energies are in MeV. For the numerical evaluations, 
we take a representative supernova at the Galactic center 
region at d = 10 kpc, with a total energy output of 3 x 
10 53 erg equipartitioned in all neutrino flavors, i.e., £ a = 
5 x 10 52 erg for each of the 6 flavors. Further, we choose 
(E) to be 12 MeV for v e , 15 McV for D e , and 18 MeV 
for the 4 other flavors represented by v x , respectively 
We show, in Sec. IIV1 how the proposed measurement is 
highly sensitive to SN emission parameters. 



B. Detection cross section 



The yield from elastic scattering on protons, i.e. 



p — > v + p 



(any flavor) 



(4) 



is comparable to that of inverse-beta reactions (p e + p — > 
e + + n). The total cross section is about a factor of four 
smaller, but this neutral current channel couples to all 
active flavors of v and v, as opposed to only v e . The 
crucial difference is that the elastic scattering events are 
at low quenched energies T' < 2 MeV, and one needs a 
significantly low threshold to detect these events [29| . 

The differential cross section da/dT for a neutrino of 
energy E to produce a proton recoil of kinetic energy T, 
to zeroth order in E/m p , is given by j29l. |36| 



da 
dT 



G%m n \ ( _ m p T 
~ 2E 2 
4.83 x 10~ 42 cm 2 



rtipT 
~2E 2 



MeV 



1 + 466 



T 
E 2 



(5) 



where T and E are in MeV and we have used m p = 938 
McV, c v = 0.04, and c a = 1.27/2 [11]. The recoil 
kinetic energy is minimum, T = 0, for a grazing collision, 
and maximum, T max , when the neutrino momentum is 
reversed. The cross section rises linearly by a factor 
of ~ 2 over this allowed range of recoil energies T = 
(0 — T max ), i.e., higher recoil energies are preferred. Note 
that the recoil direction cannot be measured, due to 
the isotropic emission of scintillation light, so one can 
determine the energy of the neutrino using the recoil 
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TABLE I: Detector properties, i.e., number of free proton targets (N p ), Birks constant (fcs), and energy resolution (AT'/T'), 
and signal yields above 0.2 MeV for the large scintillator detectors considered here. Note that the Birks formula for quenching 
in the KamLAND detector includes the quadratic correction, while others do not. 



Detector 


Mass 


Chemical composition 


N p 


k B 


AT'/T' 


Signal Yield 




[kton] 


(rounded to nearest %) 


[10 31 ] 


[cm/MeV] 


(T' in MeV) 


(T' >0.2 MeV) 


Borcxino 


0.278 


C9H12 


1.7 


0.010 


4.5%/ VT 7 


27 


KamLAND 


0.697 


Ci 2 H 26 (80%v/v)+C 9 H 12 (20%v/v) 


5.9 


0.0100 


6.9%/VT 7 


66 


SNO+ 


0.800 


C6H5C12H25 


5.9 


0.0073 


5.0%/VT 7 


111 



energy only on a statistical basis. The cross section 
for antincutrinos is slightly different, but in practice 
this difference is negligible at SN neutrino energies [29| . 
These differences due to weak magnetism corrections 
almost cancel between neutrinos and antineutrinos, and 
become only important when average neutrino energies 
are > 30 MeV |29|, so we ignore them in our analysis. 



C. Detector response 

We estimate the number of free proton targets using 
the fiducial mass and composition as 



N A AlY^ 



Wi fi 



6.02 x 10 32 • iV 



Wifi 



)>>, 



(6) 



where the fiducial detector mass M (in ktons) is com- 
posed of different components i with weight fractions Wi, 
molecular weights rrii (in a.m.u), each contributing /, free 
protons per molecule, and Na is the Avogadro number. 

A proton with recoil energy T loses energy by repeat- 
edly colliding with electrons in the detector material. 
The rate of energy loss is pre dicted to be approximately 
oc 1/T by Bethe theory [37|. The energy loss rate is 
about (10 3 -10 2 ) MeV/cm for the considered range of 
recoils with T = (1 — 5) MeV, which is much more 
than the typical 2 MeV/cm for a relativistic electron 
in a carbon target. Also, protons have velocities in 
the range j3 = (0.03-0.07), which at the lower end are 
comparable to atomic electron velocities and the Bethe 
approximation is no longer valid (see [38| for details). 
These subtleties are accounted for by using accurate 
numerical tables for (dT/dx) taken from the PSTAR 
tables at http://physics.nist.gov. Similar data is 
also available at http : //sr im . org. 

The energy loss on hydrogen targets is significantly 
larger (almost a factor of 2) than that on carbon, so we 
account for the composition of each detector and add the 
(dT/dx) for protons on carbon and hydrogen targets in 
the ratio of their weights in the detector. 



While all the recoil energy is deposited in the detec- 
tor, only part of it leads to scintillation light. This 
"quenching" is an important effect because the proton 
is slow. The quenching function T'(T) maps a recoil 
kinetic energy of the proton T to an electron-equivalent 
quenched kinetic energy T' as 



T'(T) 



dT 



l + k B (dT/dx) 



(7) 



where k B is Birks constant |39| . This parametrization for 
the quenching function is approximate. An additional 
term in the denominator, C (dT/dx ) , is sometimes 
included to achieve a better fit to data [40| . For a typical 
Birks constant of 0.01 cm/MeV, there is no quenching 
when (dT/dx) < 100 MeV/cm. Whereas for (dT/dx) > 
100 MeV/cm, the light yield is quenched. 

Although the relationship between T and T' is non- 
linear, it is one-to-one and can be calibrated accurately. 
For typical parameter values noted above, the highest 
recoil energies T ~ 5 MeV are quenched by a factor 
T/T' ~ 2, whereas at T ~ 1 MeV, the quenching 
factor is ~ 5. Clearly, quenching affects the recoil 
spectrum and the number of signal events above a given 
threshold. See BFV [II] for details. Of the detectors we 
consider, KamLAND has published its quenching factor 
including the quadratic correction, while Borexino and 
SNO+ quote their Birks constant, or equivalent. 

The energy resolution of the detector depends on the 
number of detected photoelectrons per unit energy, i.e., 
drip e /dT'. Assuming dn pe /dT' to be almost constant 
in the relevant regime, one gets (dn pe / dT')T' hits at 
energy T' . In the limit of large number of photoelectrons, 
their Poisson fluctuations are well approximated by a 
Gaussian, which leads to an energy resolution 



AT'/T' = \/^j(dn pe /dT> 



T 



(8) 



The yield is (dn pe /dT') w few hundred detected pho- 
toelectrons/MeV, leading to a resolution AT'/T' better 
than 10% in the relevant energy range above 0.2 MeV. 
We simulate the energy resolution by smearing the signal 
locally at each energy with a Gaussian of width given by 
the energy resolution. 
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D. Backgrounds 

The backgrounds at a scintillator detector can arise 
from either steady detector backgrounds, or a variety of 
other charged/neutral current signals due to the SN itself. 

Cosmic ray induced backgrounds are very small over 
the relevant times [29j], and the most important steady 
backgrounds come from radioactivities in the scintillator 
and surroundings. Of these, the most obvious is f3- 
decays of 14 C that produce a high rate of electrons below 
0.2 McV. This is a common background at all carbon- 
based scintillators, and sets a threshold, below which 
the signal is almost completely background dominated. 
Pulse shape discrimination may be used to reject this 
background and lower the threshold, probing lower 
energy neutrinos and greatly enhancing the yield, but 
we do not assume that here. 

Above 0.2 MeV, most background events are due to a- 
decays of 210 Po in the energy range T = (0.2-0.5) MeV. 
In their common quest to detect solar neutrinos, all 
the detectors we consider have purified their scintillator 
to similar levels, and the 210 Po rates are similar and 
manageably small. This background has been measured, 
including its spectral shape, and can be subtracted 
statistically or by pulse shape discrimination. 

Charged current signals from SN neutrinos that could 
be important backgrounds at these energies are cither 
small, or can be tagged and subtracted [29j. Other 
neutral current channels, i.e., reactions on 12 C, have a 
high energy threshold [29[ and are disfavored, relative 
to BFV [20| , in the light of low average neutrino energies 
that are preferred by state-of-the-art simulations [13, El- 
Further study of the inelastic neutral current channels are 
needed [2{|. Elastic scattering events on carbon nuclei 
are heavily quenched and and unobservable at present 
detectors [29[ . Finally, elastic scattering on electrons has 
a small rate and produces larger recoil energies. 

Consequently, we can use an almost universal descrip- 
tion of the detector backgrounds for the SN burst signal - 
a sharp threshold at T' ~ 0.2 MeV, and a well understood 
reducible background above threshold. This was hoped 
for in BFV, and it is remarkable that it has been achieved 
in working experiments. 




Quenched Kinetic Energy T ' [MeV] 

FIG. 1: Galactic SN neutrino-proton elastic scattering event 
spectrum at Borexino. The dashed vertical line at 0.2 MeV 
shows the threshold used. 



A. Signals 

With these detectors, the expected neutrino signals in 
the v-p elastic scattering channel are shown in Fig. Q] 
for Borexino, in Fig. [5] for KamLAND, and in Fig. [3] for 
SNO+. We find that all three detectors can detect a 
significant number of events, i.e., ~ 100/kton, with a 
modest number of background events, above 0.2 MeV. 
We have checked that our results agree in detail with 
BFV [2t| when identical inputs are used. Note that most 
of the signal in the v — p elastic scattering channel comes 
from the hotter v x flavors. This is partly due to their 
average energies being higher than those of v e and v e , 
and partly because there are 4 flavors that contribute 
to the signal. Also note, as an interesting aside, that 
the total yield with zero threshold is about ~ 500/kton, 
similar to the v e yield from inverse beta decays at water 
Chcrcnkov detectors. In the remainder of this section, 
we justify our choices for the detector parameters. 



Borexino 



IV. DETECTED SIGNALS 

We consider the large scintillator detectors, i.e., Borex- 
ino and KamLAND, which are already available now, 
and SNO+, which should be operational shortly. For 
each of these detectors, we need to know three relevant 
quantities: (i) the number of free proton targets N p , 
(ii) Birks constant fcsj and (Hi) the energy resolution 
AT'/T'. The values of these detector parameters are 
summarized in Tabic HI 



A SN signal will be localized in time and allow use of 

0. 278 kton as the fiducial mass [U [42[ instead of the 
smaller volume used for solar neutrinos. The detector 
material is pseudocumene (C9H12) [IH, with a specific 
density 0.875. The number of free proton targets is thus 
1.7 x 10 31 . We could not find a direct published reference 
to the Birks constant for the Borexino scintillator. 
However, we determine ks = 0.010 cm/ McV using the 
available quenching data on protons at Borexino 44| . 

1. e., protons of energy 8.3 MeV and 4.6 MeV get 
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FIG. 2: Galactic SN neutrino-proton elastic scattering event 
spectrum at KamLAND. 



FIG. 3: Galactic SN neutrino-proton elastic scattering event 
spectrum at SNO+. 



quenched to 4.1 MeV and 1.86 MeV, respectively and 
a particles of 5.4 MeV are quenched by a factor of 
13. Note that this is lower than the ks found for 
electrons in Borexino [45| because low energy electrons 
are quenched more than protons or ions of the same 
kinetic energy [46[. About 500 photoelectrons/MeV are 
detected at Borexino [IJ , which sets its energy resolution 
to be AT'/T' = 4.5%/x/T 7 (with V in MeV). Finally, 
the backgrounds at Borexino at energies below 0.2 MeV, 
due to 14 C, make events at those energies unusable. 
Above that threshold, the background is negligible. It 
is comprised of the entire event yield in ~ 10 s without 
a SN burst (47[, and there are ~ 2 background events, 
mostly from the 210 Po peak. 



KamLAND 

The Kamland fiducial volume is assumed to be 
0.697 kton [48| for a SN burst, corresponding to the 
inner 5.5 m sphere of the detector. The detector 
material is a mixture of 80% (by volume) of dode- 
cane (C12H26) with 20% pseudocumene (C9H12). The 
number of free proton targets is thus 5.9 x 10 31 . 
From the observed quenching of protons with recoil 
energies (1 — 10) MeV, the quenching factor was re- 
ported Hi to be 1 k B = (0.0100 ± 0.0002) cm/McV and 
C = (2.73 ± 0.08) x 10~ 5 (cm/MeV) 2 . The energy res- 



1 Ref. |49H reports ks in units of g/cm 2 /McV which we have 
converted to units of cm/Me V using known specific densities for 
dodecane (0.750) and pseudocumene (0.875). 



olution at KamLAND is determined by a photoelectron 
yield of 210/MeV, giving AT'/T' = 6.9%/ \/T [50|. Note 
that these are marginally different from the values in 
Ref. [49|, for the pre- purification scintillator. Again, the 
backgrounds at KamLAND at energies below 0.2 MeV, 
due to 14 C, make events at those energies unusable. 
Above that threshold, the background is comprised of 
the entire event yield in the absence of a SN burst 51 1, 
wherein one expects ~ 31 background events 50J, 
mostly from the 210 Po peak. The background rates 
are known very well, and the expected fluctuations in 
the background are relatively small (~ 5 events), so we 
expect that these events can be statistically subtracted. 



SNO+ 

The SNO+ fiducial volume is taken to be 0.8 kton. 
The detector material is linear alkyl benzene 
(C6H5C„H2„+i). The alkyl group typically has a 
size n = (10 — 16), and we assume n = 12 for 
definiteness and a specific density of 0.86. The number 
of free proton targets is thus 5.9 x 10 31 . In recent 
laboratory tests of the detector material, the Birks 
constant for the scintillator in SNO+ is reported to be 
about 0.0073 cm/MeV [52j . The energy resolution at 
SNO+ is expected to be 5.0% at 1 MeV and 3.5% at 
3.4 MeV respectively [53j], from which we estimate an 
energy resolution AT'/T' = 5.0%/Vf 7 . These values 
could change in the full detector. A measurement of 
backgrounds at SNO+ is not available yet, but SNO+ 
has solar neutrino physics goals that are similar to 
Borexino and KamLAND, so we expect that the region 
above 0.2 MeV to have similar backgrounds. 
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FIG. 4: Sensitivity (defined as fractional change in event yield 
relative to benchmark SN values) to SN emission parameters 
with KamLAND as the detector. 



FIG. 5: Sensitivity (defined as fractional change in event yield 
relative to benchmark values for the KamLAND detector) to 
detector parameters for the benchmark SN emission. 



B. Sensitivity and Robustness 

Wc now investigate the expected signal dependence on 
different choices of SN neutrino fiuence parameters and 
detector characteristics. To quantify this sensitivity, we 
define S as the fractional change in the event yield for 
non-benchmark values for the SN or detector parameters. 

In Fig. |4[ we show the signal sensitivity to the average 
energy of v x , overall shape, and total energy. For these 
calculations we have assumed the detector properties to 
be that of KamLAND. Using the SN model described 
in Sec. I11II as the benchmark, we show how the signal 
varies for (i) a v x average energy (E„ x ) =21 MeV or 
16 MeV, (ii) a Maxwell-Boltzmann spectrum with the 
same average energy (E Vx ) = 18 MeV, and (Hi) a higher 
or lower v x total energy, E Vx = 1.3 or 0.7 times the 
benchmark value. We can see that the expected signal 
depends strongly on the average energy or the spectral 
shape. This is because the signal is mainly from the 
"tail" of the distribution, which is exponentially sensitive 
to (E). The dependence on £ Vx , which sets the overall 
normalization, is linear, as expected. The yield varies 
by a factor of few in the range of SN emission values 
commonly predicted from SN theory [Tol . ITll ] . 

In Fig. [5j we show the dependence on relevant detector 
characteristics. We choose (i) a higher and lower 
quenching, i.e., ks = 0.0102 cm/Me V with C = 
2.81 x 10- 5 (cm/MeV) 2 , and k B = 0.0098 cm/MeV 
with C = 2.65 x 10~ 5 (cm/MeV) 2 , respectively, i.e.. the 
reported range of uncertainty in ks for the KamLAND 
scintillator [49j| . (ii) a worsening of the KamLAND 
energy resolution to AT' /T' = 10% /VT 7 . We find 
that these uncertainties leads to very small effects on 



the predicted event yield. The signal varies by only a 
few percent in the range of uncertainties of the detector 
parameters. 

We checked the dependence on energy loss rate (not 
shown) by using a 1/T fit to the Bcthe theory prediction, 
instead of accurate numerical tables. Such a naive choice 
of (dT/dx) leads to an erroneously larger yield (~ 50% 
more). It is therefore important that accurate values for 
(dT/dx) be used. 

It is thus clear that in the range of detector un- 
certainties the signal changes by few %, while the 
statistical uncertainties with ^100 detected events are 
about 10%. The expected signal variation between 
allowed SN emission models is much larger, with changes 
of up to a factor of two. Therefore this signal can be 
expected to identify specific emission scenarios. 



V. EXTRACTING THE NEUTRINO FLUENCE 

We have shown that the quenched recoil spectrum of 
protons due to neutrinos from a SN burst can be reliably 
measured at scintillator detectors. The measurement is 
almost completely background-free above 0.2 MeV, the 
detector properties arc well-calibrated, and the energy 
resolution doesn't play a crucial role. The task at hand 
is to invert this signal into a measurement of the neutrino 
fiuence as a function of energy. 

One could choose a parametrization for the spectrum 
and fit for those parameters using data. Two obviously 
interesting quantities are the average energy (E Vx ) and 
the total energy E Vx . In BFV, fitting for those parame- 
ters, with about ~ 100 signal events, led to an expected 
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precision of about 1 / ~ 10%, though the parameters 
are correlated [29j ] . 

We show here that it is possible to take a more 
general approach. One can extract the spectrum non- 
paramctrically, and the extracted spectrum has high 
fidelity to the true spectrum, even for realistic choices 
of detector characteristics. 

The problem of reconstructing the neutrino spectrum 
from the proton recoil spectrum would be quite simple 
if there was a known one-to-one relationship between 
the recoil energy and the neutrino energy. In that case, 
as for inverse-beta events, one could have simply scaled 
the observed number of events by the cross section and 
detector size, and translated the observed energies T' to 
neutrino energies using their known relationship. 

This does not work for elastic scattering events. The 
differential cross section is broad and slightly forward 
peaked [2{|, but the recoil energy T is uniquely related 
to the minimum neutrino energy E m [ n = *J m p T/2 that 
is capable of producing that recoil. Furthermore, T is 
related in a unique way to the quenched recoil T 1 . So, 
having measured T' , we can uniquely determine T using 
the known quenching function T'(T), and then T leads to 
its corresponding E m i n . This allows us to use the recoil 
data as a function of quenched recoil energies to extract 
the neutrino fluence as a function of neutrino energy E. 



A. Inversion recipe 



Suppose the quenched recoil data is in a range of 
energies which we divide into -/Vbin bins. We denote 
the value of T" at the midpoint of each bin by T[ and 
the width of the bin by AT/, with i going from 1 to 
TVbin- Each T" is uniquely related to some recoil energy 
T = T(T') using the inverse of the quenching function. 
Each of these recoil energies T in turn are related to a 



minimum neutrino energy Ej = yJm p T(Tj)/2. 

The data is simply the number of observed events iVj 



in the i th bin 



N,=AT![^L 
dT' 



(9) 



Using the expression for dN/dT', this can be written as 

(10) 



N bin 

N i= J2 KijFj , 

3=1 



where 



Fa 



(dF/dE) Ej AE J , 

dT 



(11) 



Ku = 



N AT' ( — i (S2.\ 





for i < j 
for i > j 



because only neutrinos with energy more than */ m p T/2 
are able to produce a proton recoil energy T. We 
illustrate the form of Eq. (fTOj) explicitly for iVbin = 3, 



(12) 



It is easy to see how the above equation generalizes to 
any number of bins. 

Now the extraction of the neutrino spectrum is simply 
an inversion of this matrix K. That is, we write 



( N1 \ 




/ K u K 12 K 13 \ 




fFl \ 


N 2 




K 22 K 23 




F 2 






{ K 33 ) 




Kf 3 J 



N,, 



(13) 



and find the neutrino fluence at the detector dF/dE = 
Fj/AEj in the energy bin around E = Ej. The inversion 
can be done trivially using back-substition, as the matrix 
Kij is upper triangular. That is, we can start with 
determining the fluence in the highest energy bin and 
proceed to lower ones. Again, with A/bin = 3 as an 
example, wc have 



F 3 = N 3 /K 33 , 

F 2 = (N 2 - F 3 K 23 ) /K 22 , 

Fx = (N 1 -F 2 K 12 -F 3 K 13 )/K 11 



(14) 



Note the upper-triangular form of the matrix . This is 



It is again easy to see how the above formula generalizes 
to any number of bins. 

The procedure is unchanged even in the presence 
of known backgrounds Bi in each bin, except that, 
after the Fi are determined, one would subtract the 
expected backgrounds. This subtraction procedure works 
to an accuracy w \fBlj(Fi + Bi), i.e., as long as the 
backgrounds are small. 

There is freedom in the choice of number and width 
of bins. The optimal number of bins depends on a 
compromise between finer sampling or lower noise in each 
bin, subject to the constraint that bins are wider than 
the energy resolution. Excessively narrow bins have too 
few events and lead to spurious features from noisy data, 
whereas excessively wide bins lead to a breakdown of 
Eq. (jlip . which uses a linear interpolation within each 
bin. We find that choosing the bins such that the events 
are almost equally distributed between them, leads to an 
accurate reconstruction. 

Finally, a short remark about the stability of our 
prescription. Formally, this inversion problem is known 
as a Volterra integral equation of the first kind. Volterra 
equations are relatively stable, because the kernel 
is upper triangular. There is, however, literature on 
how to numerically stabilize the matrix inversions, if the 
kernel is even mildly singular [54| . As a safeguard against 
unphysical artifacts in the reconstructed spectrum, due 
to noisy data and possible instability of the inversion 
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FIG. 7: Reconstructed neutrino fluence from the mock 
data for KamLAND (upper panel) and residuals from the 
benchmark model (lower panel). Note the unequally sized 
bins for the data. The fiuences are to be compared to Eq. $J$, 
accounting for {E Vw ) and £ Vx . 



procedure, the inversion procedure may need to be 
regulated. The regulated solution is more stable, and 
smoothed over expected statistical fluctuations. For our 
purposes, we find that we don't need a regulator when 
we bin the data as above. 



B. Numerical example 

To demonstrate the viability of this procedure we 
generate mock data at the KamLAND detector, using 
our benchmark SN model. The data (66 signal events) 
falls in the range V = (0.2-5.9) MeV, which we divide 
into 8 bins with comparable numbers of events in each 
bin. The lower end of the observable range in T" is set 
by the large 14 C background below 0.2 MeV. The upper 
end is where the neutrino fluxes become negligible. The 
chosen bin widths are significantly larger than the energy 
resolution and thus no significant correlation is expected 
between bins. 

These data, without statistical noise, are shown in 
Fig. [5J with the number of events in each bin mentioned 
alongside. Expected one-sigma errors in the reconstruc- 
tion due to Poisson fluctuations arc plotted as error bars. 
We have already subtracted the modest backgrounds. 
This figure has a logarithmic scale on the abscissa, and 
unequally-sized bins. On a linear scale, all the bins would 
have comparable areas. 

We apply our inversion procedure to the mock data 
and recover the total fluence dF/dE. We assume that 
the v e contribution will be known accurately using the 
data from the inverse-beta channel, while v e events are 
negligible. Thus, we subtract the subleading contribution 
expected from i> e , and divide by 4 to find the v x fluence 
in each flavor - , v r , , and v T . 

The range of available recoil energies maps to neutrino 
energies in the range E = (25-73) MeV. The lowest 
available energy is set by backgrounds, and is close to the 
typical average neutrino energies expected. The higher 
end goes well into the tail. The reconstructed v x fluence 
in this energy range is shown in the top panel of Fig. [7] 
using the points for discrete energy bins, with their one- 
sigma fluctuations shown with error bars. Although we 
have shown here a reconstruction for mock data without 
statistical fluctuations, we have separately checked that 
a reconstruction for noisy data typically remains within 
the shown error bars. In the inset, we show the region of 
the neutrino spectrum that is being probed. Note that 
the peak energies are somewhat lower than the average 
energies, due to the asymmetric spectral shapes. We also 
plot some alternate fiuences - one with an average energy 
energy (E Vx ) =21 MeV, one with a Maxwell-Boltzmann 
spectrum with (E Ux ) = 18 MeV, and one with a higher 
total energy £„ x being 1.3 times the benchmark value, to 
test the reconstruction. 

The reconstructed spectrum shows that only the higher 
part of the v x spectrum is probed, because v x energies 
lower than ~ 25 MeV are swamped by background. This 
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is clearly a drawback of this channel. Notwithstanding 
this, it must be emphasized that for determining non- 
thermal features, the tail is in fact the most interesting 
region of the spectrum [55[ . The spectrum near the peak 
does not depend as strongly on average energy, and would 
not have significantly improved the ability to distinguish 
between spectral parameters. Also, signatures of flavor 
mixing in the signals, due to charged current events, 
appear at the higher energy tails where the flavor spectra 
of v e and v e are most different from that of v x . Thus it 
is quite important to be sensitive to the tail of the SN 
neutrino energy spectrum. 

The main advantage of this more general reconstruc- 
tion is that we obtain a direct measurement of the v x 
flavor spectra in the higher energy range. This data 
is obtained without any fitting or parametrization of 
the primary neutrino spectra, and can be used as an 
independent measurement to compare with the high 
energy v e or D e spectra measured in charged current 
interactions, to see if mixing occurs. Also, this provides 
a model independent comparison with SN simulations. 
Any departure for thermal spectra, or any unexpected 
features, can be seen directly. 

Note that knowledge of the distance to the supernova 
is needed to correctly normalize the neutrino fluence. 
However, it is not crucial for determining the average 
energy, or for discovering nonthermal features. On the 
other hand, for an accurate measurement of the total 
energy, or for comparing with other flavors, the distance 
must be determined by other means. 

In the lower panel, we plot the residuals from the true 
spectrum. With about 10 events in each bin, we expect 
1/VlO ~ 30% uncertainties in each bin. Our bins are 
chosen to be wide enough to be almost uncorrclatcd. On 
the other hand, when a parametric expression for the 
spectrum is being tested, the statistical power in different 
bins can combine and lead to smaller uncertainties of 
about l/\/l00 ~ 10%. These are the capabilities of 
a detector like KamLAND. If one combines the results 
from multiple detectors, e.g., KamLAND, Borexino, and 
SNO+, the uncertainties are 30% smaller. 

The events in the tail have a lot of statistical power to 
distinguish between models. In particular, the final data 
point that is far out in the tail seems to be very powerful, 
but it is also subject to larger systematic errors. This is 
simply because the highest energy event decides the bin- 
width of the last bin, and large Poisson fluctuations there 
can lead to an wrong estimate of the average flux in that 
bin. With actual data, the binning will need to make 
optimal use of the available signal. We find that, even 
after omitting the final bin, the alternate scenarios are 
disfavored with A% 2 > 3.5. 

Note that the reconstructed spectrum is not identical 
to the true spectrum. This is due to biases in our 
discretization. We evaluated the cross section in each 
bin at the midpoint. Given that the differential cross 
section has a term that varies as 1/E 2 with E, it means 
that we assumed a larger than effective cross section. On 
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the other hand, the flux itself is falling exponentially in 
E. Our prescription to calculate all quantities in a bin 
at the mid-point has a bias for the largest sized bins. 
For the SN signal, we find that the bias in our scheme 
is within statistical uncertainty, as long as the binning is 
not too crude. The bias can be reduced by considering a 
properly weighted scheme, as opposed to the mid-point 
scheme we have employed. 



VI. CAPABILITIES OF A LARGER DETECTOR 

In this section we discuss how the potential of this 
channel increases if a larger scintillation detector is 
built. To illustrate the potential, we choose LENA, a 
proposal for a (50-100) kton scintillator detector as a 
possibility (BH . \5l\ . The number of free proton targets 
is expected to be 3.3 x 10 33 , which corresponds to a 
fiducial mass of 44 kton We assume the Birks 

constant for the scintillator to be 0.010 cm/MeV, in the 
ballpark of other scintillators. The energy resolution is 
expected to be worse than smaller experiments, due to 
more absorption in scintillator of the scintillated light, 
and we assume AT' /T' = 10%/ s/T 7 . These values 
are likely to be different in the future detector, but we 
believe that our choices are conservative. Backgrounds at 
LENA are not available yet, but LENA has solar neutrino 
physics goals that are similar to Borcxino , KamLAND, 
and SNO+, so we expect that the region above 0.2 MeV 
to be similarly background free for a SN burst. With 
these assumptions, and our benchmark SN model, the 
expected event rate at LENA is shown in Fig. [8] 

We find approximately 5250 events, almost comparable 
to the yield of v e at Super-K for a similar SN burst. With 
such a large number of events, we expect this data will 
also allow for a time-dependent analysis. This may reveal 
important time dependent variations in the v x emission 
properties. These possibilities need to be investigated in 
more detail, but we do not pursue them here. 

Applying our reconstruction procedure to this data 
from a LENA-typc experiment will lead to extraordinary 
results. One could opt for finer binning in energy, as 
long as allowed by energy resolution, or simply much 
smaller statistical uncertainties in each bin. We choose 
the former option and, in Fig. |9l show a reconstruction 
with 16 bins. With about 70 times larger statistics than 
KamLAND, we get a precision of < 5% in each bin. 
Thus we expect spectral and luminosity parameters to 
be determined to a precision of few %. Although biases 
begin to become important, they remain manageable. 
The statistical power is significantly better than what 
can be achieved at present detectors, and would give 
strong constraints on acceptable SN neutrino spectra, 
once systematic uncertainties are under control. 



VII. WHAT CAN WE LEARN? 

The most important message is that the v-p elastic 
scattering channel is potentially accessible at already 
available scintillator detectors, and the triggers required 
to gather the relevant data must be put in place. Missing 
out on this crucial piece of data would be a huge loss to 
SN neutrino phenomenology. 

The primary neutrino physics result from this mea- 
surement would be a direct measurement of the v x 
spectrum. This information complements the flavor 
information available at water Cherenkov and liquid 
Argon detectors. This would allow for easy comparison 
between these detectors and one would be able to identify 
the flavor exchanges in both the "disappearance" and 
"appearance" channel, providing a complete picture. 

A revealing aspect would be the relationship between 
the observed u ei D e , and v x spectra. First, one would 
check if there are non-thermal features in the observed 
spectra. If so, one would ask if the non-thermal features 
appear at identical energies for different flavors. An 
answer in the affirmative would reveal that flavor con- 
versions are at the heart of this observed non-thermality. 
This would reveal the pattern of flavor exchanges over the 
observed energy range, shedding light on the unknown 
neutrino parameters, i.e., sign(Am^ tm ) and sin #13. 
Additionally, this allows one to test the equipartition 
hypothesis and compare with SN simulations. 

This data would for the first time allow us to empiri- 
cally test the claim that almost all the SN energy goes in 
neutrinos. Only when we have detected all flavors can we 
find the energy output in neutrinos and compare that to 
the binding energy released by the star. Determination 
of the binding energy will also be a useful diagnostic for 
the proto-neutron star mass and radius, which relates to 
the neutron star mass and radius [59146 1| . The measured 
v x spectrum is also a probe of nucleosynthesis [62T - l66j . 

Any oscillation into sterile neutrinos would be most 
readily observable using this neutral current channel. 
This will set stringent bounds on their masses and 
mixing, through constraints on total energy loss and a 
flavor composition in different detection channels. 

It is important that all three experiments - Borexino, 
KamLAND, and SNO+, be running and actively looking 
for these events. The most obvious reason is to increase 
statistics. However, having more than one detection 
also allows for useful cross-calibration to rule out any 
unexpected backgrounds. Additionally, a multi-detector 
search makes it unlikely that this important signal would 
be missed owing to down-time for a specific experiment. 

Even the smallest of the considered experiments, 
Borexino, is expected to see ~ 27 events, and is capable 
of measuring the total energy at about 20%. Combining 
all three available detectors gives ~ 204 events above 
threshold, with ^ 10% precision on total energy. 
Average energies that are significantly larger ( & 10%) 
than our benchmark value would be easily distinguished 
or constrained. 
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The proposed large liquid scintillator detector LENA 
could be extremely useful for this channel. At present, 
its detector specifications are not yet certain. However, 
with conservative estimates on the specifications, we find 
that the prospects are promising. One can expect ~ 5250 
events, and few % level measurements of SN v x emission 
parameters. This is almost comparable to what Supcr- 
K can do for v e . Clearly, such a significant detection in 
multiple flavors will allow for meaningful comparison of 
the data from these two detectors. 



VIII. SUMMARY AND OUTLOOK 

We have presented an updated calculation of the SN 
neutrino signal due to elastic scattering on protons at 
scintillator detectors. Using more realistic assumptions 
on SN emission and detector properties than were used 
in BFV, we find that the signal is observable at available 
detectors. Additionally, we have demonstrated a simple 
procedure to numerically reconstruct the v x spectrum 
from the data, which will allow for new analyses and 
probe different aspects of SN physics and astrophysics. 

As a note for the future, we must remark that detec- 
tors with lower quenching will lead to more promising 



results. Another avenue for drastic improvement is if the 
threshold can be lowered below the 14 C background. This 
may be possible through pulse shape discrimination [44| . 
That allows us to extend our range to lower neutrino 
energies, increases the overall yield drastically, and may 
allow reconstructing the SN spectral peak. Together, one 
could hope for significantly better results than we have 
outlined for already achieved experiments. 

In conclusion, we exhort the experimentalists working 
in the Borexino, KamLAND and SNO+ collaborations, 
and on proposed detectors like LENA, to seriously 
consider the physics impact of this channel and be 
adequately prepared to acquire this data from a future 
Galactic SN signal. 
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